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Abstract 

Mass transfer from an evolved donor star to its binary companion is a standard feature of stellar 
evolution in binaries. In wide binaries, the companion star captures some of the mass ejected in a wind 
by the primary star. The captured material forms an accretion disk. Here, we study the evolution of 
wind-accretion disks, using a numerical approach which allows us to follow the long term evolution. 
For a broad range of initial conditions, we derive the radial density and temperature profiles of the 
disk. In most cases, wind-accretion leads to long-lived stable disks over the lifetime of the AGB donor 
star. The disks have masses of a few times 10"^ — 10~^ -Af©, with surface density and temperature 
profiles that follow broken power-laws. The total mass in the disk scales approximately linearly with 
the viscosity parameter used. Roughly 50% to 80% of the mass falling into the disk accretes onto 
the central star; the rest flows out through the outer edge of the disk into the stellar wind of the 
primary. For systems with large accretion rates, the secondary accretes as much as 0.1 M©. When the 
secondary is a white dwarf, accretion naturally leads to nova and supernova eruptions. For all types 
of secondary star, the surface density and temperature profiles of massive disks resemble structures 
observed in protoplanetary disks, suggesting that coordinated observational programs might improve 
our understanding of uncertain disk physics. 



1. INTRODUCTION 

Mass transfer from an evolved donor star to a com- 
panion star is a common outcome of stellar evolution in 
binary systems. In close binaries with periods of 1-2 yr 
or less, the donor can fill its tidal surface and transfer 
mass rapidly into a massive disk surrounding the com- 
panion (e.g.. iKenvon et al. 1982; Bat h fc Pringlelll982bl : 
iKenvon et al.lll991l : iSiviero et al. .20091 ). Although these 
disks are lui ninous, the systems are short-lived and ra rely 
observable (jKenvon fc Webbink|[l98llKenvonlll986D . 

Disks formed via wind accretion may be more com- 
mon. Chemically peculiar stars, including Barium stars 
and CH stars, are probabl y pollut ed by material accreted 
from a binary companion (IMcClur c & Woodsworth 1990; 
ILuck fc Bondlliggit IHan et all 11995; .Busso et al..,2001, ). 
White dwarfs (WDs) accreting material through disks 
filled by the winds of their evolved companions are 
often identified a s symbiotic stars or very slow (also 
symbiotic) novae (IKenvon fc Truran"1983: Ke nvoniri986l : 
IWebbink et all Il987t ISokoloski et al.. .2006.) . In both 
cases, long-lived disks are possible when the donor star 
does not fill its tidal lobe. 

Theoretical studies of the physical properties of wind- 
accretion disks provide an important framework for inter- 
preting observations of these systems. The replenishment 
of newly accreted material and the formation of a disk 
may also have important implications for the formation 
and evolution of planetary systems. A newly formed disk 
may interact with pre-existing planets, leading to the re- 
growth and/or migration of these planets, or p ossibly 
the fo rmation new planets and/or debris disks (jPeretsI 
I2010allbh . 

Several previous studies explored the formation of 
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wind-accretion disks through SPH and AMR simulation s 
(jMastrodemos fc Morrislll99 8': 'de Val -Borro et al.ll2009l) . 
Here we use a different simplified analytical/numerical 
approach typically used for the study of disks in shorter 
period binaries and in protostellar/protoplanetary sys- 
tems. This approach allows us to follow the long term 
evolution of the disks, which is more difficult in hydro 
simulations, and to compare the derived physical struc- 
ture with protoplanetary disks. We focus on the accre- 
tion from low mass companions in wide binaries (3 — 100 
AU) , which evolve through wind-accretion. 

We begin with a description of the basic conditions in 
wind-accreting binaries, followed by a detailed explana- 
tion of the methods we use to follow their evolution. We 
then describe results of the wind- accretion disk structure 
followed by a discussion and summary. 

2. WIND ACCRETION IN WIDE BINARIES 

To understand the conditions required for the forma- 
tion of ci rcumstellar disks from win d accreted material, 
we follow lSoker fc RappaportI (|2000f ) . For a disk to form, 
ja , the specific angular momentum of the accreted mate- 
rial, must exceed J2 = (GAf2^2)^^^, the specific angular 
momentum of a particle in a Keplerian orbit at the equa- 
tor of the accreting star of radius R2 ■ For accretion from 
a wind, the net specific angular momentum of the ma- 
terial entering the Bondi-Hoylc accretion radius, Ra^ is 
JBH = 0.5(27r/P)i?2 (Wang 1981), where P is the orbital 
period. The actual accreted specific angular momentum 
for high Mach number flows is ja = iiJbh, where 77 ~ 0.1 
and 77 - ^0.3 for isotherm al and adiabatic flows, respec- 
tively (jLivio et al.lll986[ ). For a binary composed of a 
main sequence (MS) accretor with mass M2 and a mass- 
losing asymptotic giant branch (AGB) star with mass 
Afi, requiring ja > J2 leads to a simple condition for 
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where R2 is the radius of the accreting star, a is the semi- 
major axis of the binary, and Vr is the relative velocity 
of the wind and the accretor. 

Adopting i?2 '^ Rq , disks can form in binaries with a < 
10-100 AU for Ml « 1-10 Mq and Vr ^ 5 - 20 km s-^. 
This range in a spans t he peak of the distributio n 
of binary star separations (iDuauennov fc MavoH ITQQII ) . 
Thus, many binaries contain wind-accretion disks during 
the late stages of their evolution. 

Red giant branch (RGB) and AGB stars lose mass in 
low velocity win ds at typical rates of M^, ~ 10~^ — 10~^ 
Mq yr-i (e.g., 'Taylor & Seaquist' 1981 jKem^ [TOSl 
lOlofsson et al. 2002; Ramstcdt et al. 200(|. To estimate 
the fraction of this material accreted by the MS compan- 
ion, we adopt the Bondi-Hoyle rate for an isotropic wind 
from the mass-losing star. 
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where Ra is the Bondi-Hoyle accretion radius 
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In this expression, v'^ — v^ 4- v^ is the velocity of wind 
material near the secondary, Vq is the orbital velocity of 
the secondary, and Cs is the sound speed of the wind. 
For typical velocities of ~ 10-20 km s~^, Ma/M^ ranges 
from ~ 20% for close binaries with a « 3 AU, to < 0.1% 
for wider binaries with a « 100 AU. Thus, maximum 
accretion rates can briefly reach roughly 10~^ Mq yr~^ 
for close binaries with AGB star companions. 

Despite the relatively low efficiency, wind accretion can 
add considerable mass to the secondary. While on the 
RGB and AGB, the primary star loses most of its initial 
mass. For primaries with initial masses 1-7 Mq, close 
companions can accrete as much as 1 Mq from wind 
accretion. 

Detailed AMR simulations of wind- accretion disks 
yield results consistent with these simple estimates 
(jde Val-Borro et aD 120091 ). When the star serves as a 
sink term in the simulations, circumstellar material forms 
a fairly stable disk with a total mass of ^ 10~^ of the stel- 
lar mass out to radii of roughly 5 AU. Typically, ^ 5% of 
the mass lost from the AGB star goes through the disk, 
consistent with the simple Bondi-Hoyle estimates. 

To develop a model for the structure of the disk, we as- 
sume the primary loses mass at a constant rate, M^ ■ For 
wide binaries, wind material is ejected far from the accre- 
tor at approximately parallel trajectories. Gravitational 
focusing close to the accretor then sets the cross-section 
for material hitting the disk. Instead of the geometrical 
cross-section, the fraction of the wind captured by the 
disk up to distance r is fj^ (< r) oc r. For an outer 
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The differential accretion rate per radial bin in the disk, 
is therefore 

Macc{K)dR = KCl^ . (5) 

The source function, the instantaneous surface density 
profile of material accreted into the disk is 
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As wind material falls into the disk, it dissipates kinetic 
energy. Together with irradiation from the central star 
and viscous heating, this energy contributes to the heat- 
ing of the disk. The amount of energy input into an 
annulus at radius r is the gravitational binding energy. 
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Although the disk does not accrete material beyond 
Ra, viscous processes can expand the disk beyond Ra- 
However, tidal forces from the binary companion limit 
the outer radius. Typically, this outer radius is roughly 
90% to 95% of the Roche limit, 02 = ria, where 
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and q = A/2 /Mi (lEggletonlll9"83[ ). 

For wind mass loss with Mi > M2 and low accretion 
efficiency, binaries expand. We approximate this expan- 
sion assuming adiabat ic mass loss from the primary star 
(|Hadiidemetrioulll963D 
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where Mi^mit are the initial masses of the binary compo- 
nents, Qinit is the initial separation and 



M,{t) = A/,, 



M^{t)dt, 
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for i = 1,2. Assuming a constant mass loss rate from the 
primary, Mi(i) = Mi,,„,t - M^t and M2(i) = M2,„,t + 
C^M^^t, where Q is the accretion efficiency. For a system 
losing mass in a wind {C, <C 1), a{t) > amit- Thus, the 
binary separation expands with time. For most mass 
loss rates and separations of interest, the disk reaches a 
steady-state on timescales shorter than the lifetime of the 
AGB star (see results in Table 1). Thus, we can safely 
adopt a simple prescription for the outer radius of the 
disk 

Rout = 0.2 (11) 

and assume that a is constant in time. Considering a 
broad range of a and M^ allows us to derive conditions 
in the disk for a broad range of plausible binaries (see 
Table 1). 



The inner boundary of the disk, Ri- 
radius of the accreting star 

Rin = R2', 



is hmited by the Solving this equation for Tyy requires 
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we adopt an inner radius of 1.5i?2; the precise value of 
the inner radius has little impact on the overall structure 
of the wind-filled accretion disk. 

Equipped with the necessary initial and boundary con- 
ditions, we now discuss a model for the evolution of the 
disk. 

3. DISK EVOLUTION 

3.1. Numerical calculation 

For a disk with surface density E and viscosity v, con- 
servation of angular momentum and energy leads to a 
non-l inear diffusion equation for the time evolution of E 
(e.g.,Eindcn-Bell & Pringlc 1974; Pringlc 1981), 
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The first term is the change in E from viscous evolu- 
tion; the second term is the change in E from other 
processes, incl uding mass loss from photoevaporation 
(e.g. [A lexande r et al.l I2Q06D or planet formation (e.g., 
[Alexander &: Armitagd I2009D and mass gain from wind 
material falling into the disk (given by Eq. [5]). Our 
approach neglects the torque of the disk by the wind, 
which is small for a disk in a binary system. For the 
wide binaries we consider (30 and 100 AU separation), 
the wind momentum can compress the outer edges of the 
disk but does not affect its bulk structure. The viscos- 
ity is i^ = aCsH, where Cs is the sound speed, H is the 
vertical scale height of the disk, and a is the viscosity 
parameter. The sound speed is c^ = jkTd/ ^imn, where 
7 is the ratio of specific heats, k is Boltzmann's con- 
stant, Td is the midplane temperature of the disk, /j, is 
the mean molecular weight, and ruH is the mass of a hy- 
drogen atom. The scale height of the disk is H = Csfl^^, 
where fl — y/GM^^/R^ is the angular velocity. 

To solve this equation numerically, we assume that the 
midplane temperature is the sum of the energy generated 
by viscous (T^) dissipation, the energy from irradiation 
(T/), and the energy of inf ailing material from the wind 

T^^n+Tf + T^. (14) 

The viscous temperature is 



n^ 



27kvY?Q'^ 



64cr 



(15) 



where k is the o pacity and a i s the Stefan-Boltzmann 
const ant (e.g. iRuden fc LinI 119861 : iRuden fc PoUackl 
[T99lh . With V = acln-^ and is = (27a/64cr) - 
{jk/ ^ttih) '^ nflT,^, the viscous temperature is 



T^ - t2Td. 



(16) 

The energy from infalling material follows from Eq. [T) 
Balancing the kinetic energy of infall with the energy 
radiated by the disk yields 
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If the disk is vertically isothermal, the irradiation tem- 
perature is Tf{R) = (ei2){R^/RfT^, where i?* and T* 
are the r adius and effective temper ature of the central 
star and (jChiang fc Goldreichiri997[ ) 
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Thus, the irradiation temperature is 
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Follow ing [Chiang fc Youdml (l2010f l and lHueso fc Guillotl 
(l2005h . we set dlnH/d\nR = 9/7. With H = Csil'^ 
we set to = {2T^/3Tr){R^/Rf and ii = [RjRf - 
[IRVty^ ~ {^k/fj.mHy/^ - T^. The irradiation tem- 
perature is then 



Tf^to + tiTy\ 



(21) 

Viscous disks are not vertical ly isothermal 

(|Ruden fc P ollack '19911 : iD'Alessio et"all ll998^: however 
this approach yields a reasonable approximation to the 
actual disk structure. 

Because Ty and T/ are functions of the midplane tem- 
perature, and Tw is independent of it, we solve equation 
(14]) with a Newton-Raphson technique. Using equations 
( TB|) . (PT|) . and[T51 we re- write equation p^ as 



f{Td) = TJ- (to + tiT]/2 



Adopting an initial Td 
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^2' "^ or Trf « t\' ' , the derivative 
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and yields a converged Td to a part in 10® in 2-3 itera- 
tions. 

In the inner disk, the temperature is often hot enough 
to vaporize dust grains. To account for the change in 
opacity, we follow IChamberj (|2009[ ) and 



assume 
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with n = -14 in reg. i ons with T^ > Te = 1380 K 
(jRuden fc Pollacklll99lt lStepinskilll998l) . For simplicity, 
we assume k = kq when Td < T^. 

To solve for the time evolution of E, we use an explicit 
technique with N annuli on a grid extending from Xm to 
Xout where x = 2 R^^"^ (jBath fc Pringl c 1981, 1982a), and 
Rin and Rout are inner and outer bo undaries determined 
by Eq s. [TT] To verify this code, iBromlev fc Kenvoiil 
(|2011h compare numerical solutions with analytic results 
from [Chambers. (,2009.V 



3.2. Potential caveats 

Our study explores wind- accretion disks using simple 
numerical models. These models provide us with the 
ability to explore the long term evolution of disks with 
a wide range of initial conditions. Because these disks 
reach an approximate steady-state, we can test whether 
the solutions scale as expected with the binary separation 
and the mass loss rate of the primary star However, the 
approach makes many simplifying assumptions, and does 
not consider many physical processes known to occur in 
these systems. Because the physics we do not include is 
often poorly understood or weakly constrained, it seems 
prudent to begin with a simple physical structure and 
add additional physics as needed to understand a broader 
range of phenomena. The following paragraphs discuss 
useful physics to consider for future studies. 

Radiative heating of the disk by the companion: 
Our calculations do not include radiative heating of the 
disk by the primary star (mass donor). If the primary has 
a radius i?i « 0.5 AU and a luminosity Li, a flared disk 
surrounding the secondary with outer radius Rout and 
vertical scale height hout ~ 0.05 R„ut intercepts 0.4% 
to 1% of Li for a = 3-100 AU (jKenvon fc HartmannI 
119870 . With ii « 500-1000 L©, this contribution to 
disk heating is comparable to irradiation by the central 
MS star. In wide binaries with a w 10-100 AU, AGB 
stars can have i?i « 1 AU. The disk then intercepts 0.2% 
(a = 10 AU) to 0.02% (a = 100 AU) of the radiation from 
the primary. For typical AGB Li ^ 10"* L0 and a « 10 
AU, this component is somewhat larger than irradiation 
from the central star. In all configurations, the outer 
rim of the disk intercepts nearly all of this radiation. 
If the outer disk radiates efficiently, this extra heating 
should not impact the structure of material in the in- 
ner disk. To test this assumption, we performed several 
test simulations with a fixed outer boundary tempera- 
ture, aT^^t = firrLi/AnRouthout- This extra heating 
reduces the disk surface density by a factor of ~ 2 near 
the outer boundary but changes the surface density at 
smaller radii, R < 0.95Rout, by less than 1%. Thus, ne- 
glecting this component has little impact on our results. 
If the outer disk radiates inefficiently, the outer disk will 
expand and transport thermal energy radially inward. If 
the scale height of this puffed up disk becomes compara- 
ble to the accretion radius Ra it could significantly affect 
the accretion process, and would require a different ap- 
proach. 

Disk viscosity: The processes underlying the origin 
of the viscosity in accretion disks are still not understood. 
Here, we follow most studies and adopt a simple prescrip- 
tion using the a viscosity parameter. In particular, we 
do not directly consider here any interactions of mag- 
netic fields and their evolution, and we use a constant 
a in each model. Our general results are shown our for 
a choice of a = 0.01. Table 1 also provides additional 
data for the cases of a = 0.001 and a = 0.1, enabling us 
to consider the overall dependence of disk structure and 
evolution on the viscosity parameter. 

Binary orbital expansion: As briefly discussed 
above, we do not change the binary separation and mass 
in our calculation. In most cases the wind-accretion disk 
evolves to a steady state in a relatively short time scale. 
Our approach then provides a good approximation to re- 



ality. When evolution in the binary separation occurs 
more rapidly than the time-scale to reach steady-state 
conditions in the disk (e.g., when the time-scale to reach 
steady-state is > 10^ yr), our approximation is less ac- 
curate. Our results should then be taken with more cau- 
tion. In general, both the evolution of the binary and 
the equilibrium disk mass evolve slowly enough to treat 
the time evolution with a sequence of steady disk models 
in binaries with increasing semi-major axis. 

Accretion model and disk instabilities: We use a 
simplified ID model for the accretion infall of wind mate- 
rial into the disk. Although this model captures many of 
the main properties of the material infall, reality is prob- 
ably more complex. Aside from our failure to treat the 
complex interaction between the wind and the outer disk, 
our simplified 1.5D model does not allow for variations 
in the vertical structure which might lea d to thermal in- 
stabi lities and dwarf nova outbursts (e.g.. lCannizzo et al.l 
I2OI2I and references therein) or for wind s driven by radi- 
ation pressure from the inner disk (e.g.. lNoebauer et al.l 
120101 ). In particular, the inner region where the tem- 
peratures rise above 1000-2000 K, (roughly where dust 
grai ns evaporate) might be come thermally unstable (see 
e.g. [Alexander et al.l 1201 1[ ). Though the location of the 



disk instability regime depends on the opacity model, the 
mass loss rate of the AGB primary, and other input pa- 
rameters. Fig. [3]shows that temperatures as high as 1000 
K are typically achieved only for a small region, extend- 
ing up to an AU for the models with the highest accre- 
tion rate and/or the smallest binary separations. In most 
models, these relatively high temperatures are achieved 
only over small regions in the inner 0.1-0.3 AU of the 
disk. In models with the smallest binary separation (3 
AU) , the outer disk rad ius is com parable to disks studied 
in other contexts fe.g.. [Alexander et al. 2011) . In wider 
binaries with separations of 10-100 AU, the outer disk 
extends well beyond regions likely capable of sustaining a 
thermal instability cycle. Thus, potentially unstable re- 
gions can extend over a large fraction of the disk in com- 
pact binaries (a w 3 AU) but impact a fairly small frac- 
tion of the disk in wider binaries (a > 10-100 AU). We 
therefore conclude that models exploring small separa- 
tions, particularly those with the highest accretion rates 
(e.g. models 1 and 5 in table 1 below), could be highly 
susceptible to disk instabilities in a significant fraction 
of the inner disk. We caution that our steady state so- 
lutions might not well represent these cases. The steady 
state solution is likely to well represent likely outcomes 
in wider binaries and in compact binaries with low accre- 
tion rates. In all cases, likely unstable regions are close 
to the accreting star and could could lead to outbursts 
with distinct observational outcomes. In Alexander et 
al. (2011), the optical brightness changes little during 
the outbursts; the bolometric luminosity changes by a 
relatively small factor of a few for the largest infall rates, 
though these are larger than those we study. 

4. RESULTS 

To develop an understanding of the long term behavior 
of wind-fed disks in wide binaries, we consider a grid of 
models in (a, Mwind) space (see Tables 1 and 2 for model 
parameters and results). For simplicity, we adopt q = 
3, M2 = 1 Mq (MS secondary) or M2' = 0.6 Mq (WD 
secondary), and assume constant separation and mass 



loss rate with time. For each (a, M^ind) pair, we set Vw 
= 10 kni s~^ and Cg — 20 km s~^, derive Ra from equation 
([3|) and Rout from equation (fTTj) . The disk inner cut-off 
is i?in = 1.4_R2, with i?2 = 1 i?o for the main sequence 
stars and R2 — 0.01 Rq for the white dwarf stars. 
For an initial disk mass Af ^ — 10~^ M^i, the initial 
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surface density distribution is E = T,oR^^e~^^^°^^ with 
So = Mdfi/'i'nRRout- Using our numerical solution to 
the diffusion equation, we continuously add material to 
the disk inside Ra and evolve the surface density in time 
until (and if) the disk reaches a steady state at time 
Tsteady, Or up to 1 — 2 Myr, corresponding to the max- 
imum lifetimes of AGB stars. As long as the steady- 
state disk mass Msteady is much larger than the initial 
disk mass, numerical tests show that the initial surface 
density distribution has little impact on the final steady- 
state surface density distribution. 

We consider the system to reach a steady state when 
the change in disk mass per unit time is smaller than 
10~^ of the input rate. In most cases, disks reach steady 
state on timescales much shorter than the AGB lifetimes 
(see Fig. [T]). Several models with the largest binary sep- 
arations do not properly fulfill this steady state condition 
by the end of the simulation, but the actual change in the 
disk structure at late times is small. For these models 
we show the disk properties at the end of the simulation 
run time, Tmn- In these cases, Trun is already compa- 
rable or longer than the typical AGB lifetime. Thus, 
wind accretion disk cannot achieve a steady state con- 
figuration. Nevertheless the disk structure in these cases 
hardly changes over most of the AGB lifetime. 

Tables 1-2 summarize the main results of the calcu- 
lations. The detailed steady state profile and evolution 
of the disks and the accretion rate (on to the star) are 
shown in Figs. [T]|31 Although we focus on main sequence 
star accretors in the Figures, results for WD accretors 
are similar. Table 2 shows results for white dwarfs with 
negligible luminosity. If the accreting star is a hot white 
dwarf with L = 100 L©, UV photons heat up disk re- 
gions close to the central star but impact the outer disk 
very weakly. Thus, the time evolution of the accretion 
rate and disk mass (in Fig. 1) and the radial profiles of 
cumulative mass, surface density, and disk temperature 
(Figs. 2-5) are fairly representative of WD accretors and 
MS accretors. 

For all calculations, the disk mass and the accretion 
rate onto the star increase roughly linearly with time and 
then reach an approximate steady state (Figure [1]). The 
timescale to reach steady-state increases roughly linearly 
with increasing binary separation. Although the relation 
is shallower than a linear one, binaries with smaller mass 
infall rates also take longer to reach steady state than bi- 
naries with large mass infall rates. The steady state disk 
masses range from ~ 10"'^ M© for binaries with massive 
winds and small separations to ~ 10^^ M© for bina- 
ries with low mass loss rates and large separations. For 
models with identical mass loss rates from the evolved 
companion, wider binaries have smaller infall rates into 
the disk and smaller disk masses. When binaries with 
different a have similar infall rates, wider binaries have 
larger and more massive disks. 

For most main sequence star models, wind-fed accre- 
tion rates produce negligible luminosities. With typical 
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Fig. 1. — Disk mass profile and evolution. Top: The evolution of 
the total mass in the accretion disk. Disks in binaries with larger 
separations are wider and accommodate more mass, as they are 
truncated at further distances from the star. Bottom: Evolution 
of the accretion rate onto the star. Both upper figures show how 
the disk evolves into a steady state configuration over timescales 
which are typically much shorter than the mass loss time scales. 
The lines correspond to the models described in Table 1. Lines of 
the same color correspond to similar mass loss from the companion 
(IQ-'^, lO-*^, IQ-^ and 10"* M© yr-^ from top to bottom). Lines 
of the same type correspond to similar binary separation (3 AU, 
dash-dotted; 10 AU, dotted; 30 AU, solid; 100 AU, dashed). 

steady accretion rates smaller than 10~^ M© yr~^, the 
accretion luminosity is comparable to or smaller than the 
stellar luminosity of 1 L©. Thus, the accretion disk and 
boundary layer are invisible (Kenyon & Webbink 1984). 
For close binaries with AGB-type primary stars, accre- 
tion rates of 10~® M© yr^^ produce modest accretion lu- 
minosities of ^ 10 \jQ. Although optical radiation from 
the disk and boundary layer are not detectable, ultravi- 
olet (U y) radiation from the boun dary layer is probably 
visible (jKenvon fc WebbinklHosl . 

Models of wind-fed white dwarf stars are more in- 
teresting. In these systems, accretion rates exceeding 
10~^ M© yr~^ yield bright UV sources which c an ionize 
the wind from the primary (jKenvon fc Webbi nk 1984). 
This ionized wi nd produ ces bright optical and UV emis- 
sion lines (Kenvon fc Web bink 1984) and a luminous ra- 



TABLE 1 

MS SECONDARY MODELS 



#i 


^,- 


^2^ 


„4 


«a 




-/o.. 


-2... 


,jO.B 9 


mI» 




^steady 


^^2„ 1 






















10-^ 


10-^ 


10-1 


10-^ 


10-^ 


10-1 


10-^ 


10-^ 


10-1 


1 


3.0 


1.0 


3.0 


1.5 


1.4 


-5.0 


-5.74 


-0.05 


-6.23 


-2.20 


-3.12 


-4.01 


4.51 


3.51 


2.57 


5.00 


4.00 


3.00 


2 


3.0 


1.0 


10.0 


2.5 


4.8 


-5.0 


-6.55 


-0.85 


-6.82 


-2.12 


-2.96 


-3.77 


5. OS 


4.28 


3.51 


5.30 


4.30 


3.70 


3 


3.0 


1.0 


30.0 


3.1 


14.3 


-5.0 


-7.46 


-1.74 


-7.60 


-2.14 


-3.00 


-3.89 


6.11 


5.28 


4.51 


6.30 


5.30 


4.70 


4 


3.0 


1.0 


100.0 


3.4 


47.6 


-5.0 


-8.49 


-2.80 


-8.59 


-2.53 


-3.49 


-4.47 


6.97 


6.30 


5.00 


6.97 


6.30 


5.30 


5 


3.0 


1.0 


3.0 


1.5 


1.4 


-6.0 


-6.74 


-1.05 


-7.23 


-3.00 


-3.85 


-4.70 


4.54 


3.67 


2.84 


5.00 


4.00 


3.00 


6 


3.0 


1.0 


10.0 


2.5 


4.8 


-6.0 


-7.55 


-1.85 


-7.82 


-2.77 


-3.62 


-4.49 


5.51 


4.66 


3.83 


6.00 


5.00 


4.00 


7 


3.0 


1.0 


30.0 


3.1 


14.3 


-6.0 


-8.46 


-2.74 


-8.60 


-2.89 


-3.82 


-4.80 


6.51 


5.51 


4.52 


6.60 


5.60 


4.60 


S 


3.0 


1.0 


100.0 


3.4 


47.6 


-6.0 


-9.49 


-3.80 


-9.59 


-3.47 


-4.46 


-5.46 


7.00 


6.00 


5.04 


7.23 


6.00 


5.30 


9 


3.0 


1.0 


3.0 


1.5 


1.4 


-7.0 


-7.74 


-2.05 


-8.23 


-3.68 


-4.52 


-5.40 


5.51 


4.51 


3.51 


5.70 


4.70 


3.70 


10 


3.0 


1.0 


10.0 


2.5 


4.8 


-7.0 


-8.55 


-2.85 


-8.82 


-3.47 


-4.39 


-5.34 


5.83 


4.96 


4.04 


6.00 


5.00 


4.30 


11 


3.0 


1.0 


30.0 


3.1 


14.3 


-7.0 


-9.46 


-3.74 


-9.60 


-3.80 


-4.77 


-5.77 


6.52 


5.54 


4.54 


6.70 


5.70 


4.70 


12 


3.0 


1.0 


100.0 


3.4 


47.6 


-7.0 


-10.49 


-4.80 


-10.59 


-4.46 


-5.46 


-6.46 


7.04 


6.04 


5.00 


7.30 


6.30 


5.30 


13 


3.0 


1.0 


3.0 


1.5 


1.4 


-8.0 


-8.74 


-3.05 


-9.23 


-4.37 


-5.28 


-6.24 


5.23 


4.36 


3.43 


5.48 


4.48 


3.48 


14 


3.0 


1.0 


10.0 


2.5 


4.8 


-8.0 


-9.55 


-3.85 


-9.82 


-4.34 


-5.32 


-6.32 


6.04 


5.04 


4.04 


6.30 


5.30 


4.30 


15 


3.0 


1.0 


30.0 


3.1 


14.3 


-8.0 


-10.46 


-4.74 


-10.60 


-4.77 


-5.77 


-6.77 


6.54 


5.54 


4.51 


7.00 


6.00 


5.00 


16 


3.0 


1.0 


100.0 


3.4 


47.6 


-8.0 


-11.49 


-5.80 


-11.59 


-5.46 


-6.46 


-7.46 


7.04 


6.00 


5.04 


7.30 


6.30 


5.30 



Columns: (1) Model number; (2) Primary mass (Mq) ; (3) Secondary mass (Mq); (4) SMA separation (AU); (5) Bondi-Hoyle accretion 
radius (AU) ; (6) Outer boundary of the disk (AU); (7) Mass loss rate from primary {Mq yr" 1) ; (8) Accretion rate into disk (Mq yr" i) ; 
(9) Total mass that goes through the disk in 0.5 Myrs (Mq); (10) Accretion rate onto the star (Mq yr" l) ; (11) Disk mass at steady state 
(Mq); (12) Time to achieve steady state (yrs); (13) Simulation run time (yrs). 
Columns 7-13 are shown in logarithmic scale. 



TABLE 2 
WD SECONDARY MODELS 



#' 


mi^ 


■J 

1712'' 


a* 


Ra^ 


Rout 


Mlo.s^ 


Md..fe« 


j.rO.S 9 
^''-'Hi.sk 


M,„io 


Md,,feii 


12 

"^steady 


T 13 

'run 


1 


1.8 


0.6 


3.0 


1.2 


1.4 


-5.0 


-6.06 


-0.37 


-6.47 


-3.29 


4.00 


4.48 


2 


1.8 


0.6 


10.0 


1.7 


4.8 


-5.0 


-6.96 


-1.26 


-7.18 


-3.24 


4.51 


4.70 


3 


1.8 


0.6 


30.0 


2.0 


14.3 


-5.0 


-7.89 


-2.19 


-8.02 


-3.25 


5.51 


5.70 


4 


1.8 


0.6 


100.0 


2.1 


47.6 


-5.0 


-8.96 


-3.24 


-9.00 


-3.51 


6.30 


6.30 


5 


1.8 


0.6 


3.0 


1.2 


1.4 


-6.0 


-7.06 


-1.37 


-7.47 


-4.00 


3.75 


4.00 


6 


1.8 


0.6 


10.0 


1.7 


4.8 


-6.0 


-7.96 


-2.26 


-8.18 


-3.85 


4.80 


5.00 


7 


1.8 


0.6 


30.0 


2.0 


14.3 


-6.0 


-8.89 


-3.19 


-9.02 


-3.92 


5.88 


6.00 


8 


1.8 


0.6 


100.0 


2.1 


47.6 


-6.0 


-9.96 


-4.24 


-10.01 


-4.39 


6.30 


6.30 


9 


1.8 


0.6 


3.0 


1.2 


1.4 


-7.0 


-8.06 


-2.37 


-8.47 


-4.64 


4.51 


4.70 


10 


1.8 


0.6 


10.0 


1.7 


4.8 


-7.0 


-8.96 


-3.26 


-9.18 


-4.48 


5.51 


5.60 


11 


1.8 


0.6 


30.0 


2.0 


14.3 


-7.0 


-9.89 


-4.19 


-10.02 


-4.70 


6.00 


6.00 


12 


1.8 


0.6 


100.0 


2.1 


47.6 


-7.0 


-10.96 


-5.24 


-11.01 


-5.33 


6.30 


6.30 


13 


1.8 


0.6 


3.0 


1.2 


1.4 


-8.0 


-9.06 


-3.37 


-9.47 


-5.28 


4.51 


5.00 


14 


1.8 


0.6 


10.0 


1.7 


4.8 


-8.0 


-9.96 


-4.26 


-10.18 


-5.21 


5.54 


5.60 


15 


1.8 


0.6 


30.0 


2.0 


14.3 


-8.0 


-10.89 


-5.19 


-11.02 


-5.60 


6.00 


6.00 


16 


1.8 


0.6 


100.0 


2.1 


47.6 


-8.0 


-11.96 


-6.24 


-12.01 


-6.29 


6.30 


6.30 



Columns: (1) Model number; (2) Primary mass (Mq) ; (3) Secondary mass (Mq); (4) SMA separation (AU); (5) Bondi-Hoyle accretion 
radius (AU) ; (6) Outer boundary of the disk (AU); (7) Mass loss rate from primary {Mq yr" l) ; (8) Accretion rate into disk (Mq yr" i) ; 
(9) Total mass that goes through the disk in 0.5 Myrs (Mq); (10) Accretion rate onto the star (Mq yr— l) ; (11) Disk mass at steady state 
(Mq); (12) Time to achieve steady state (yrs); (13) Simulation run time (yrs). 
Columns 7-13 are shown in logarithmic scale. 



dio sources at cm wavelengths (jTavlor fc Seaauistlll984D . 
These system s would be classified as symbiotic stars 
(|Kenvod ll98ffl. 

Despite their much lower accretion rates, the central 
star accretes a larger fraction of infalling material in wide 
binaries. For a Ki 3-100 AU, the Bondi-Hoyle accre- 
tion radius Ra exceeds Rst the stagnation radius where 
the flow yelocity in the disk changes sign. Thus, some 
fraction of infalling material flows out through the disk. 
The rest accretes onto the central star. In these simu- 
lations, the ratio of accreted to ejected material ranges 
from roughly 55% for close binaries (a = 3 AU) to ~80% 
for wide binaries (a = 100 AU). The accretion flow in the 
accretion disk is bimodal. Generally, material at separa- 
tions somewhat shorter than Ra flows into the star; the 
rest of the material outflows out into the extended disk 
up to Rout and is eyentually ejected from the disk. 

Accreted material typically has a negligible impact on 
the central MS star. For typical lifetimes of 10^ yr or 



less, the central star accretes less than 1% of its initial 
mass. In the closest binaries with the highest mass loss 
rates, the accretor may accrete from 10% to 90% of its 
initial mass over the lifetime of the companion (e.g. see 
the overall accreted mass assuming steady state, over 
0.5 Myr, M^dk'^ Table 1). Ahhough our models do not 
account for the mass increase of the secondary in these 
cases, the disk achieves steady-state long before the MS 
star accretes a significant amount of mass. Thus, our 
steady-state disks still correctly represent the likely phys- 
ical state of the system. 

Accreted material has a more significant impact on 
accreting white dwarfs. Aside from their visibility 
as symbiotic stars, wind-fed disks surrounding white 
dwarfs can pr oduce a variety of eruptive phenomena 
(jKenvonl I1986D . At low accretion rates, cold white 
dwarfs under go degenerate shell flashe s and become clas- 
sical novae (jKenyon fc TruranI Il983[ ). White dwarfs 
accreting at larger rates yield non-degenerate flashes 
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Fig. 2. — Radial profile of wind accretion disks at steady state. 
Top: The radial density profiles of wind accretion disks under var- 
ious conditions. Also shown is a shaded region corresponding to a 
range of protoplanetary disk models which only differ in their over- 
all normalization corresponding to the range of total disk masses 
as inferred from observations (^ 2 X 10~^ — 10~^ M0; mean of 
2 X 10-3 M,7>: [Andrews fc W illiams 2007). Bottom: The cumula- 
tive mass radial profile of wind accretion disks. Lines correspond 
to the same models in Fig. 1 (described in Table 1). 

as in the symbiotic stars V1 016 Cyg and HM Sge 
(jMikola iewska & Kenyon 1 9921 ). In some systems, the 
rates are sufficient to allow the white dwarf to evolve 
close to th e Chandrasekhar limi t and become a type la 
supernova (jStarrfield et al.|[2005D . 

Steady-state surface density profiles generally follow 
expectation. Over the entire disk, the surface density 
approximately scales with radius as S oc r"", with n « 
0.8-1.5. The normalization of the surface density scales 
almost linearly with the accretion rate into the disk. 
Changes in slope typically occur in the inner disk, when 
the disk becomes hot enough to evaporate dust grains, 
and in the outer disk, when the energy generated by irra- 
diation or i nfall domin ates the energy generated by vis- 
cosity (e.g.. IChambersl[20 09'. and references therein). 

We briefly explored the effect of different viscosity pa- 
rameters on the overall structure and mass of the disks 
(see table 1). In general, the steady state disk mass is 
roughly inversely proportional to a, and scales linearly 
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Fig. 3. — Radial temperature profile of wind accretion disks (in- 
set shows a zoom up region; note linear distance-scale). The radial 
temperature profiles of wind accretion disks under various condi- 
tions are shown; the lines correspond to similar models as described 
in Fig. 1. Note that regions in which temperatures rise beyond 
approximately 1000 K might be susceptible to thermal disk insta- 
bility, not modeled here. 

with a~^ . These results hold as long as the disk is cold. 
As the disk gets hotter, more of the disk has an opacity 
dependent temperature which slightly changes the rela- 
tion, but the overall behavior is close to linear. Beside 
the different overall scaling, the disk spatial structure is 
very similar. 

Steady-state temperature profiles are generally shal- 
lower than steady-stea dy accretion d isks fed by a lobe- 
filling companion (Lvn den-Bell fc Prin glc 1974). In most 
steady disks, the effective temperature of the disk scales 
with radius as T oc r"™, with m « 0.75. For wind-fed 
disks in very wide binaries at very low accretion rates, the 
steady-state slope approaches the standard m « 0.75. As 
the binary separation contracts and the wind accretion 
rate grows, the temperature profile becomes more shal- 
low. For the closest binaries, our results suggest m k, 
0.4. The range in the slope is small, from m k, 0.4 in 
close binaries with large accretion rates to m « 0.7 in 
wide binaries with small accretion rates. Thus, the tem- 
perature profile scales very weakly with the separation 
or the mass loss rate of the primary star. 

5. DISCUSSION 

The results from the calculations lead to several broad 
conclusions. Physical properties of wind-fed accretion 
disks depend on the orbital separation and the evolu- 
tionary state (mass loss rate) of the primary star. When 
the secondary is a main sequence star, it is observable 
only when the binary is close and the primary is an AGE 
star. Otherwise, wind-fed accretion has few, if any, di- 
rect observational consequences. When the secondary is 
a white dwarf, the binary is almost always observable as 
a symbiotic star. 

Although optical/UV spectroscopic observat ions are 
suffic ient to detect luminous symbiotic stars ([Kenvonl 
Il986f ). direct imaging can often reveal the binary com- 
panion. In Mira B, lUE spectroscopic observations 
first identified emission from a hot white dwarf accret- 



ing m aterial from the Mira wind (jReimers fc Cassatellal 
Il985f ). recent HST observations resolve d the system 
and detected the accretion disk directly (jlreland et alj 
l2007f l. For a derived binary separation of ~ 70 AU, 
the inferred accretion lum inosity of ^ lO~^''M0yr~^ 
(jSokoloski fc Bildstenll20ldl) is roughly consistent with 
our models. High resolution imaging of other symbiotic- 
like binaries, including R Aqr and CH Cyg, would pro- 
vide additional tests of these calculations. 

Aside from producing interesting symbiotic stars, WD 
accretors in wide binaries can accrete enough mate- 
rial to approach th e Chandrasekhar limit (Table 2; 
IStarrfieldet all [20051 IShen fc Bildstml l20?y7l and refer- 
ences therein) . If these WDs grow , they are sus ceptible 
to traditional type la supernovae (jArnettl [19961 ) and to 
sub-Chandrasek har He lium detonation supernovae (e.g. 
(iWooslev et all [198 6: Bilds ten et all l2007t iPerets et al] 
[20101 : iWaldman et al...2011. ). In most simulations, the 
growth of the white dwarf depends on the way mate- 
rial is accreted. Thus, understanding the structure of 
the disk is important for understanding whether wind- 
fed white dwarfs are good SNc progenitors. Our disk 
models provide a step along the path to understanding 
these outcomes. 

Wind-fed disks around white dwarfs provide a less 
traditional connection to other astrophysical systems. 
For separations of 3 — 100 AU and mass loss rates of 
10~* — 10^^ ^Q^' wind-fed disks have surface density 
and temperature profiles similar to those observed in 
low-mass protoplanetary disks. For comparison. Fig. [2] 
shows the range of surface density profiles inferred for an 
ensemble of protoplanetary disks in nearby star-forming 
regio ns ([Andrews & Williams 2007; Chiang & Youdin 
[20101) . The shaded region corresponds to disks with 
masses of - 2 x 10"'' - 10"^ Mq. With typical hfe- 
times of a few Myr, physical processes in protoplanetary 
disks yield planets on short timescales. Given the similar 
lifetimes of AGB mass-losing primary stars, it is plausi- 
ble that some aspects of planet forma tion occur in the 
wind- fed disks of wide binary systems (|Peretsl[2010al[bl ). 

Despite the similarity in instantaneous disk masses, 
disks in wide binaries have a distinct advantage over pro- 
toplanetary disks in producing planets. Continuous feed- 
ing by the primary star guarantees that the total mass 
available for planet formation is, in principle, larger for 
wide binaries than for protoplanetary disks. In close bi- 
naries with evolved AGB primaries, the total mass ac- 
creted by the central star is at least as large as the most 
protoplanetary disk. Thus, massive planets could grow 
in many of these disks. 

Testing this idea is challenging. The most massive 
disks in wide binaries are not as large as those in pro- 
toplanetary disks; dust emission from the AGB primary 
complicates acquiring the submm and mm observations 
required to estimate dust masses. Occasional nova erup- 



tions may also frustrate the coagulation processes that 
grow dust grains from the Mira wind into the planetes- 
imals that produce planets. However, ALMA has suffi- 
cient resolution to detect structures in the disks of wide 
binaries. Identifying structures similar to those in pro- 
toplanetary disks might allow robust comparisons be- 
tween the grain properties and the profiles of surface den- 
sity and temperature. Detection of debris disks around 
nearby white dwarfs suggests that planets (or debris from 
planets) survive the evolution from a main sequence star 
into a white dwarf (e.g., Kilic et al. 20ll). The frequency 
of wide WD companions to WD debris disks provides 
some estimate for the likelihood of this 'second genera- 
tion' planet formation scenario. 

6. SUMMARY 

Using a set of numerical calculations, we have explored 
the evolution and long term steady state structure of 
wind-accretion disks in wide evolved binaries (a w 3 — 
100 AU). These systems evolve rapidly (see Fig. [T]) and 
achieve a steady state on timescales which are typically 
much shorter than typical stellar evolution timescales of 
r^ 10^ yr. During steady state, disks have similar surface 
density and temperature profiles with total masses of a 
fewxlO"^ — 10"'^ M0. The radial density profile is a 
broken power law in the inner regions and an exponential 
decline at the outer edge where the disk is truncated by 
tidal forces from the companion. The radial temperature 
profile is also described by a broken power law but does 
not decline dramatically at the outer edge of the disk. 

Understanding the formation, structure, and evolution 
of wind- fed accretion disks is important for a wide variety 
of phenomena in evolved wide binary systems, including 
chemically peculiar stars, novae, supernovae, stellar out- 
bursts, and symbiotic binaries. In close binary systems 
with evolved AGB primaries, significant accretion from a 
wind-fed disk provides a natural mechanism for chemical 
enrichment and abundance anomalies in otherwise nor- 
mal main sequence and white dwarf stars. Over a broad 
range of separations and mass loss rates, wind-fed accre- 
tion onto white dwarfs can produce nova eruptions and, 
possibly, supernova eruptions. By providing a founda- 
tion for relating the mass loss rate of evolved red giants 
to the temperature and luminosity of the companion, our 
results also enable more detailed studies of the physical 
structure of individual binary systems. 

Our analysis begins to make a link between the disks in 
evolved binaries and protoplanetary disks. The structure 
and evolution of both types of disks depends on uncer- 
tain physics, including the disk viscosity and interactions 
between stellar photons and disk material. For most ac- 
cretion rates, large disks often drive massive winds which 
may interact with the surrounding wind from the pri- 
mary (in a binary system) or a molecular cloud (in a 
protostellar system). Observational comparisons among 
disks with similar accretion rates and sizes might improve 
our overall understanding of accretion phenomena. 
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